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Abstract 

In the recent paper |j^] the order a 4 Roc corrections to the positronium P levels 
were reconsidered. Those calculations confirm our corresponding results, except for 
the contribution due to the squared spin-orbit interaction. We present here a new 
derivation of our previous result for this last correction, this derivation being to our 
opinion both simple and convincing. 
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High-precision measurements of positronium level structure provide a unique test of 
quantum electrodynamics. Of course, as accurate theoretical calculations are necessary for 
this test. In Refs. || || we have found analytically the corrections of order a 4 Poo to the 
positronium P levels. 

As mentioned in |2|, ||, this theoretical problem for states of nonvanishing angular 
momenta is quite tractable, the main difficulty being of the "book-keeping" nature. 

One could be only pleased that this tedious calculation has been checked in recent paper 
0. Our results for almost all contributions to the correction discussed, are confirmed 
there. The only exception is the term induced by the spin-orbit interaction squared. The 
corresponding operator is 

v =£-4 {slpxp,])2 - (1) 

Here S is the total spin operator, p and p' are momenta of the initial and final states, 
respectively; q 2 = (p' — p) 2 . Our result for the correction induced by this operator is 

5E = tn 96^0 v 1 " i) [2(1S)2 + (1S) + 2S{S + (2) 

Here e n = ma 6 /n 3 , m is the electron mass, n is the principal quantum number, 1 is the 
orbital angular momentum. The result of Ref. |I[ is different: 

8E = e n JL (l - JjL) [2 (1S) 2 + (IS) + 2S(S + 1)] - e n ^S(S + 1). (3) 

The disagreement is most probably due to different treatment of the matrix element 
in the coordinate representation at r — > 0, which is the only subtle point in the whole 
problem. Therefore, after getting acquainted with Ref. [|I|, we did not confine to performing 
again the calculations in the coordinate representation in two independent ways, though 
this check confirmed our formula (0) and in particular the fact that for a P-wave both 
structures 2(1S) 2 + (IS) and 2S(S + 1) should enter with the same coefficient. We used as 
well a different approach, where the most part of calculations was done in the momentum 
representation, which allows one to get rid at all of the spurious singularities at r — > 0, or 
at large q. This is the approach we present below. 

The momentum part of the expectation value we are interested in, can be written as 

a f dp f dp' 4n 

km = 8^J (27p^ (p) tlkl e i mnPlPn J J2^f^ m ^ (P } - (4) 

Here 

^(p) = ^epF(p), W ) = \/JeVF(p') 

are the wave functions of P-states in the momentum representation, e, e' being constant 
unit vectors. The integral over p' can be presented as 

f dp' Ail 

J ~(^Kf~q^ V ' iV ' iV ' k F<kP '^ = A ( p } piP i Pk + B (p)( 6 ijPk + $ikPj + SjkPi)- (5) 
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The only term in the rhs contributing to the correction (Q), is B(p)5ijPk- To calculate it 
explicitly, let us multiply Eq. (||) by (5y — PiPj/p 2 )pk- in this way we obtain 

B(P) = ±J (^^(PP')[^ 2 " (PP'Wl *V). (6) 
As to the matrix element itself, it equals now 

= 8^ / -(^ F (P)( e P)l 5 kmP 2 -PkPm](e'p)B( P ) 

x{ei[44nAj - ^(5 jm - 5 jfc 5j m ]e^}. (7) 



With the identity 

— / dQ(en)ljli(e'n) — — [ dQ[e x n]j [e' x n],- = — e^e' (8) 

we arrive at the following expression for the energy correction: 



Cy 4-7T 

^ = -^o^ ls ) 2 + ( 1S ) + 2S ^ S + ^ < 4 ^ 2 + - p ' 2 ) 2 > • ( 9 ) 

Thus, without ever running into singular expressions or integrals, we have demonstrated 
that 2(1S) 2 + (IS) and 2S(S + 1) enter the result in a sum only. It is convenient now to go 
over in the perfectly convergent radial expectation value < 4nq 2 + ^{p 2 — p' 2 ) 2 > to the 
coordinate representation. With the obvious identity 

3 



< PP ' >- -H/t""'' 

4:71 



the correction is rewritten as 



^ = "32^ [2(1S)2 + (1S) + 25(5+1)] 



, „, , ma 4 

-6\R'(0)\ 2 + — <- 4 > 



(10) 



This is a simple matter now to reproduce our formula @. 
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